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Long-term  Global  HYCOM  Objective 

o  To  depict  and  predict  the  3D  global  ocean 
state  at  fine  resolution  (0.08°  on  the  equator, 
~7  km  at  mid-latitudes) 


Strategy 

o  Begin  with  0.72°  resolution  global  HYCOM  and 
optimize  the  KPP  mixed  layer  performance  (the 
subject  of  this  presentation) 

o  Before  the  end  of  FY03,  start  running  and 
development  of  0.24°  global  HYCOM 

o  Submitted  DoD  HPC  Challenge  proposal 
(FY04-06)  to  continue  development  at 
0.24°  and  start  0.08°  model  development 


0.72  degree  Global  Domain 


Pan-Am  Global  Grid 

o  0.72  degree  equatorial  Mercator  78S-47N 
o  Arctic  bi-polar  patch  above  47N 

*  Low  resolution  global  had  patch  at  59N 

*  Can’t  include  Hudson  Bay 

o  Double  latitudinal  resolution  near  the  equator 
o  Halve  latitudinal  resolution  in  Antarctic 

Coastline  at  50m  isobath 

o  Closed  Bering  Strait 

o  No  Sigma  (terrain-following)  vertical  coordinate 
Same  26-layers  as  0.08  degree  Atlantic 
First  Z-level  3  m,  increases  1 .1 25x  up  to  12m 


47N:  SCPY  (km) 


47N:  Grid  Aspect  Ratio  (SCPX/SCPY) 


0.72  degree  Global  Standard  Configuration 

KPP  mixed  layer 
Energy-Loan  ice  model 

Sigma-theta  (some  Sigma-2  runs,  not  discussed  here) 
Horizontal  diffusion  chosen  to  suppress  eddies 
Initialize  from  Levitus 

ECMWF  Reanalysis  monthly  mean  forcing 

o  Plus  6-hrly  wind  anomalies  from  sep94-sep95 

Annual  mean  of  10  largest  rivers  via  precip  bogus 

Strong  relaxation  to  monthly  Levitus  SSS 

o  Necessary  to  prevent  SSS  drift 

Weak  relaxation  to  monthly  Levitus  SST 

o  Parameterizes  SST’s  effect  on  longwave 
radiation 


Longwave  Radiation  and  SST 

Longwave  Radiation  is  sum  of: 

o  Upward  blackbody  longwave  radiation 

*  Qbb  =  -0.98  (5.67  x  10“8)  [Ts  +  273.16)4 

o  Downward  atmospheric  longwave  flux 

*  Highly  dependent  on  cloudiness 

*  Unknown  dependence  on  SST 
(assume  independent) 

If  longwave  was  calculated  using  a  SST  of  Tso: 

°  Qiw(Ts)  =  Qiw(TS0 )  +  Qbb(Ts)  —  Qbb(Tso) 

°  Qiw(Ts)  =  Qiw(TS0 )  +  Qbb  (Ts  —  Tso ) 
o  Q'bb  =  -0.98  (5.67  x  10"8)  4  (Ts  +  273.16)3 

Or,  approximately  (in  W/m2)\ 

O  Q'bb  =  -4.506  -  0.0554  Ts 

So,  —5.3W/m2/degC  is  a  good  longwave  correction 

o  Equivalent  to  3.5  m  reference  thickness  and 
30  day  e-folding  SST  relaxation 
o  lOx  weaker  than  typically  SSS  relaxation 

Ocean  Model  Intercomparison  Project  includes 
(Qbb(Ts)  -  Qbb(TS0))  as  a  longwave  correction 


Variable  Turbidity 


Penetrating  solar  (shortwave)  radiation  is  important 
for  accurate  SST 

HYCOM  parameterizes  turbidity  via  a  “two  band” 
version  of  Jerlov  water  types 

o  One  of  5  classes,  red  and  blue  bands 
o  Same  everywhere  in  time  and  space 

NLOM  used  “single  band”  kPAR  turbidity 

o  Photosynthetically  Available  Radiation 
o  kPAR  from  SeaWIFS  k490 
o  monthly  climatological  fields 

Single  band  approach  only  OK  for  bulk  mixed  layer 

Added  two  band  Jerlov-like  kPAR  scheme 

o  kPAR  isn’t  a  good  fit  to  Jerlov 
o  Much  better  than  single  class  everywhere 

Need  a  better  scheme: 

o  3-band  based  on  chlorophyll 
(Morel  &  Antoine,  1994)? 

o  Jerlov-like  based  on  k490? 


Annual  SeaWiFS  KPAR 


SST  Metrics  I 


Run  for  5  years  and  form  monthly  means 

o  May  not  be  long  enough 
o  Length  of  run  vs  confidence  in  results 
o  Takes  two  days  on  64  IBM  POWER4  opus 

Compare  monthly  SST  to  Reynolds  and  Smith 
(R&S)  climatology 

o  Monthly  anomalies 
o  Annual  mean  difference 
o  RMS  difference 
o  Correlation  Coefficient 
o  Skill  Score 

*  Correlation  squared  -  Unconditional  Bias 
-  Conditional  Bias 

*  Maximum  is  1,  but  minium  is  -infinity 

*  Measure  of  error  w.r.t.  seasonal  cycle 
(i.e.  w.r.t.  standard  deviation) 

*  Use  a  minimum  of  1  degC  for  standard 
deviation 

•  Still  get  poor  skill  scores  near  equator 


SST  Metrics  II 


Purpose  of  comparison  is  to  find  “good  enough” 
configuration 

o  Assume  that  “skill”  on  climatological  forcing  is 
maintained  on  interannual  forcing 

o  Is  monthly  thermal  climatological  forcing  enough? 

o  NLOM  experience  suggests  that  this  is  OK,  but 
can’t  be  certain  until  we  run  interannually  with 
HYCOM 

Targets: 

o  Annual  mean  error  <  0.5  degC 
o  Correlation  Coefficient  >  0.6 
o  Skill  Score  >  0.3 

Use  zonal  averages  to  reduce  amount  of  data 

o  Average  not  necessarily  best  statistic 

*  A  few  large  negative  skill  scores  can 
dominate  the  average 

o  Same  targets  as  for  full  field 


Simulation  History 


Expt  1 .5: 

o  Standard  KPP 
o  Jerlov  la  everywhere 
o  No  rivers 
o  No  SST  relaxation 

Expt  1 .8: 

o  Monthly  kPAR  turbidity 
o  Annual  rivers 

Expt  3.3: 

o  “Longwave”  SST  relaxation 
Expt  3.9: 

o  Modified  KPP  shear  instability 
Expt  4.2: 

o  Twin  of  3.3  (same  winds) 
o  FNMOC  thermal  forcing  (average  98-01) 
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ZONAL  SST  MEAN  ERROR,  HYCOM  VS  Reynolds&Smith 
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ZONAL  SST  METRICS  (III),  HYCOM  VS  Reynolds&Smith 
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ZONAL  SST  METRICS  (III),  HYCOM  VS  Reynolds&Smith 
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03.3  vs  R&S  SST:  Mean  Error 


20E  40 E 


04.2  vs  R&S  SST :  Mean  Error 


03.3  vs  R&S  SST:  Correlation 
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03.3  vs  R&S  SST:  Skill  Score 
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04.2  vs  R&S  SST:  Skill  Score 


Conclusions 


KPP  is  performing  well  in  HYCOM 

Western  Equatorial  Atlantic  a  major  trouble  spot 

o  Modifying  KPP  shear  instability  helps 
o  More  needs  to  be  done 

In  general,  most  of  the  SST  error  is  in  the  annual 
mean 

Not  yet  clear  how  much  is  due  to  forcing  and  how 
much  due  to  KPP 


